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Abstract. Quantum transport through an open periodic array of up to five dots is 
investigated in the presence of a magnetic field. The device spectrum exhibits clear 
features of the band structure of the corresponding one-dimensional artificial crystal 
which evolves with varying field. A significant magnetically controlled current flow 
is induced with changes up to many orders of magnitude depending on temperature 
and material parameters. Our results put forward a simple design for measuring with 
current technology the magnetic subband formation of quasi one-dimensional Bloch 
electrons. 



PACS numbers: 73.23.Ad, 73.21.La, 73.20.At, 73.23.-b 
1. Introduction 

Single quantum dots are tlie solid state analogue of an atom whereas the properties 
of coupled-dots may resemble that of molecules. Quantum transport through open 
quantum dots, being an equally intriguing as well as extensively investigated topic [1], 
continues to provide new insights in fundamental phenomena and fuels a wealth of 
nanoelectronic applications [1-8]. Arrays of coupled-dots may be considered as one- 
dimensional artificial crystals with the dot repeating unit acting as the lattice basis. If 
the coupling is strong enough the formed Bloch states yield an electronic structure that 
uncovers many similarities with the subbands of quasi one-dimensional systems with a 
much reduced reciprocal lattice vector. It is also well known that a uniform magnetic 
field applied to Bloch electrons yields magnetic subbands with an overall different 
spectrum [9, 10]. Unlike the lack of any impact in one dimension, in two dimensions these 
form the famous Hofstadter butterfly [11]. The effect of confinement in two dimensional 
ribbons has been studied in Ref. [12]. The question remains open as to what extent there 
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exists an observable magnetic effect for tfie intermediate dimensionality, as in the case 
of an array of open quantum dots. Moreover, experimental evidence in the literature 
is scarce [13, 14] and the effect of magnetic subbands is hard to isolate in the common 
setup of lateral semiconductor superlattices [15]. Hence, the prospect of measuring its 
properties in a simple fashion is quite attractive. 

In this study we consider small coupled-dot arrays that present distinct spectral 
properties regulated via an applied magnetic field B. The electron transport exhibits 
bright and dark windows reflecting an electronic structure that is reminiscent of the 
energy bands of the corresponding linear artificial crystal. This unique feature allows 
to explore the S-dependence of the subbands of the quasi one- dimensional Bloch 
electrons. With varying magnetic field, our calculations demonstrate qualitative (and 
quantitative) changes of the bright and dark transport windows in the suggested array 
structure, thus, yielding a direct signature of the magnetic subband formation in the 
magnetoconductance. 

Coupled-dot arrays may also be used as elements in magnetoresistive devices [16]. 
Such a design of either chaotic [7] or rectangular [17] quantum dots in alignment has 
been recently realized with a split-gate technique. In particular, the experiments of 
Ref. [17] showed a large magnetoresistance at a field slightly greater than the magnetic 
field Be that corresponds to an electron cyclotron radius equal to the size W of the 
dot, i.e.. Be = hkp/eW {kp is the Fermi wavcvcctor). Despite the possible role of 
Bragg rcficction in a periodic array, in these experiments it is still unclear how just 
Bragg reflection of electrons would result in the non-exponential drop of the conductance 
measured as the number of the dots increases. In fact, new conductance features such as 
enhanced reflection can be expected ior B/Bc^ 1 purely from the induced changes in the 
classical dynamics due to commensurability between the cyclotron radius and W. Here, 
by excluding any other contribution to the magnetoconductance, we expose the pure 
quantum mechanical effect of Bragg reflection of electrons propagating phase-coherently 
across the dot array. Magnetic subbands form for any B/B^ < 1, thereby, suppressing 
the conductance due to Bragg reflection at the newly formed band edges. Hence, as B 
varies, the flngcrprints of single bands are exposed via a significant magnetoresistance 
even at moderate field strengths. 

2. Formulation of the problem 

Fig. 1(a) shows the setup in discussion. We assume that square quantum dots 
of size W arc laterally confined on the surface of a semiconductor heterostructure 
by an electrostatic field which creates effective hard wall boundaries for ballistically 
propagating electrons. The coupled leads are modelled by quasi one-dimensional 
conduction band electrons freely propagating along the x - direction with a Fermi 
distribution fxiE) = [exp{ ^~^r^ ) + l]~^, jjix = Ef^cYsd/'^ being the chemical potential 
in the left (K=L) and right (K=R) lead when a bias voltage Vsd is applied. The point 
contacts bridging the dots have square geometry of dimensions Li, — D — 0.3W that 
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Figure 1. (Color online) (a) Schematic representation of the discussed open array 
of quantum dots, (b) Upper panel: field-free quantum transmission through a single- 
dot (dashed curve) and the five-dot array of (a). Lower panel: energy spectrum of 
the corresponding one-dimensional artificial crystal with lattice spacing L = W -\- Lf,. 
Note that flat energy bands do not contribute to transport since electrons acquire zero 
group velocity, (c) Same as (b) for a magnetic flux $ w 0.7(/>o piercing the unit cell. 
We recall that the integer part oi kpD/Tr indicates the number of propagating channels 
in the leads and q defines the Bloch vector of the periodic structure. 



are of the order of the Fermi wavelength Xp = 2TT/kF- Although quantitative details 
differ, our conclusions are independent of this simplest design. 

We model the electronic structure via a single-band effective mass equation of 
electrons in a magnetic field. The corresponding Hamiltonian is if = ^^27r^^ — ^(^) 
where m* is the effective mass (fixed to O.OSme unless otherwise stated). The boundary 
conditions are imposed via the confining potential U{r), which is chosen to be zero inside 
the enclosed area in Fig. 1(a) and infinite outside. The Hamiltonian is discretized on a 
two-dimensional tight-binding grid using the Peirls' substitution for the vector potential 
A and can be expressed in the second quantization form [18]: 

H{r) = er4cr + ^(V^e'^'^^cK+Ar + h.c). (1) 

r r,Ar 

Here, the cj (cr) define a set of creation (annihilation) operators on each site of the 
grid and Ar indicates the vector from the site r to its nearest neighbors. The on-site 
energy is er = 4:V with the hopping matrix element V = h'^/2m*a'^; a is our lattice mesh 
constant. The magnetic field B = Bz, which is applied in the region // of Fig. 1(a) 
and is zero in regions / and ///, is introduced via the vector potential A in the Peierls 
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phase factor; (pQ — h/e is the flux quantum. Charge transport properties are calculated 
within the Landauer scattering approach which expresses the current as follows: 



in which T{E) is the quantum transmission function for injected electrons with energy E; 
the factor two accounts for spin degeneracy. We calculate T using our parallel algorithm 
of the recursive Green's functions method [19]. As the system size increases one needs to 
invert a block-tridiagonal matrix which scales linearly with the array length. For serial 
processing this yields an additional cost that we avoid by distributing the scatterer's 
domain over several processors. 

3. Results 

The upper panels of Figs. 1(b) and 1(c) show the field-free and B 0, respectively, 
quantum transmission in the first open channel. Transport through a five-dot array 
is indicated by the solid curves. In contrast to the single-dot transmission spectrum 
- plotted in dashed line - bright and dark windows are formed in which transport is 
either allowed or suppressed. These compare well with the energy bands and gaps of 
the electronic structure of the corresponding infinite linear artificial lattice obtained 
independently by direct diagonalization; this is shown for zero and finite B in the 
lower panels of Figs. 1(b) and 1(c), respectively. Also evident in those figures is 
the prominent dependence of the band structure with respect to the magnetic fiux 
piercing the unit cell. Broad energy bands contribute electron states that are almost 
fully transmitted, whereas, narrow sections exhibit weaker transmission signals. The 
remarkable characteristic is that such a transmission spectrum is rapidly obtained for 
a quantum dot array with just a few unit cells as can be seen from the comparison 
of the upper and lower panels of Figs. 1(b) and 1(c). In practice, this facilitates the 
realization of a device at length scales comparable to the electronic phase coherence 
length at finite temperatures so that quantum features do not wash out. It may seem 
natural that only a small number of dots can reproduce the signatures of the artificial 
crystal, which is fully consistent with the fast convergence of the transmission with 
length in similar systems for B = (see for example Ref. [20]). It is less straightforward 
to account for the interplay of magnetic flux with lattice periodicity for flnite flelds but 
the evident correlation of energy and transmission spectra strongly supports the above 
physical interpretation of our numerical results. 

In Fig. 2, we plot the transmission function difference between the field-free 
structure and that at a field of strength B = O.SBc- The positive and negative 
parts refiect the newly formed magnetic subband structure of Bloch electrons in 
the corresponding one- dimensional artificial crystal which cause the bright and dark 
transport windows to occur at different spectral positions. As discussed later, for a 
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Figure 2. Magnetically controlled flow demonstrated via the profile of the difference 
of the quantum transmission for the field-free and B = 0.3Bc cases. 

Table 1. SI units at kpD/TT = 1.5 assuming m* = O.Qbme- 
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given geometry and Fermi energy (i.e., fixed kpO/Tr) the contrast in current flow due 
the differing transmission spectra can also be traced as a function of magnetic field to 
yield the evolution of the magnetic subbands. We note that there exists broad energy 
ranges over which bright transport windows at non-zero magnetic field overlap with 
dark areas at vanishing B, e.g., at kpD/n ~ 1.5 and kpD /n pa 1.67. This feature marks 
a mechanism for magnetically controlled current flow which can be realized up to the 
order of liquid nitrogen temperatures, as will be shown below. 

Although the system of natural units is practically convenient when estimating the 
upper limit of the conductance as determined by the number of open channels kpD/TT, 
which is independent of the number of magnetic subbands that may occur in the dot 
array, at this point it is instructive to interpret it to SI units. Assuming = 30nm with 
m* — 0.05me yields Ep — 33meV and Bf. — 1.68T. Regarding dimensions each quantum 
dot should be ~ 75nm wide and the width of the lead D 22nm at kpD/Tr ^ 1.5. 
The lattice spacing L is around lOOnm deflning a total array length of less than 500nm 
for flve coupled-dots. In a strict sense, these dimensions deflne the range of validity of 
our results regarding temperature. Apart from the thermal broadening, the temperature 
controls the scattering mechanisms determining the electronic phase coherence length. 
Since we have so far assumed that electrons are coherently propagating, the array length 
must be shorter than the latter [21]. More examples are presented in Table 1. These 
show the interplay between linear dimensions and Be- 

In Fig. 3(a), we furnish our observations with the linear- response magnetocon- 




Figure 3. (a) Linear-response magnetoconductance at various temperatures 
{kpD/'K w 1.5). Inset: magnetic field dependence of the distance between the Fermi 
energy Ep ~ 74.5mcV and the band edge Ec accounting for the resonant structure 
of the low-temperature magnetoconductance when crossing occurs at S « 0.12Bc ~ 
0.45T. (b) Typical magnetoconductance modified by the presence of weak boundary 
disorder with relative strength ai = 0.027(W) and a2 = 0.061 (W). 



ductance curve at various temperatures. An overall increase of the conductance with 
magnetic field strength is clearly observed. A remarkable feature, central to this work, 
is the fine peak-structure of the magnetoconductance dl/dV at very low temperatures 
which relates to the formation of the spectrum of Bloch electrons in a magnetic field. 
This is demonstrated in the inset of Fig. 3(a). As the band structure modifies with the 
magnetic field, the edge of a single band Ec crosses the Fermi energy Ep at B / Be ^ 0.12. 
When the distance \Ec — Ef\ vanishes a bright transport window is induced that gives 
rise to the resonant structure of dl/dV in the sub-Kelvin regime (thick fine in Fig. 3(a)). 
Due to the well-pronounced peaks one could think of using these as a probe for the mag- 
netic subband structure. At higher temperatures thermal broadening causes averaging 
over a larger part of the spectrum including many adjacent minibands and gaps. This 
increases the low-field conductance whereas simultaneously decreases the corresponding 
higher field values. 

Below we would like to discuss some aspects regarding the effects of disorder in our 
samples. Realistic devices may be experimentally prepared with electronic mean free 
paths of the order of several /^m for the regime of lower temperatures discussed herein 
[23] . These clearly allow for ballistic transport experiments. However, to account for a 
possible residual departure from the ideal periodicity of our structure we investigate two 
samples whose boundaries along the current direction are roughened with a standard 
deviation a; in this case is much larger than the disorder's characteristic length 
scale. The magnetoconductance is shown in Fig. 3(b). For a2 = 0.061(W^) the disorder- 
potential imposes a fiuctuation to the electron's energy, of ~ 2%Ep in the region of 
the square dots and ~ 30%£'j7 in the interdot bridge sections. The latter order of 
magnitude is consistent with the background disorder observed in the experimental 
setup of Ref. [17]. The presence of the disorder considered here modifies transport 
through the coupled-dot system in a two- fold way. First, the disorder within the square 
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dots' sections perturbes the position of the bands such that they may coincide with Ep, 
as it is the case in the magnetoconductance curve corresponding to the sample with 
(72. Second, the enhanced residual disorder at the bridge sections weakens the coupling 
between the dots, thereby reducing the volume of the transport windows. Effectively, 
this fact decreases the bands' width and respectively the amount of the transmitted 
current. However, the overall resonant structure of the magnetoconductance curve is 
preserved. 




B/B 



Figure 4. Linear-response magnetoconductance at T = O.IK and /cf-D/tt 3.5. 
Insets depict the band structure close to the Fermi energy at two selected magnetic 
field strengths = 0.9 (upper) and ^2 ~ 0.14 (lower) corresponding to two distinct 
conduction regimes with 1 < dl/dV < 2 and 2 < dl/dV < 3, respectively. The dashed 
vertical lines indicate the Fermi energy. 

Fig. 4 shows the conductance profile at Ep corresponding to a higher number 
of propagating modes in the leads yielding an upper limit of three. The 

magnetoconductance reveals a much richer structure. However, it is uniquely 
distinguished into regions where the maximum attainable value of the conductance 
differs depending on the number of the dot-array bands that the incoming electrons 
can populate. As illustrated using the insets in Fig. 3, with varying magnetic field more 
magnetic subbands are formed at the fixed Ep which can be further occupied by the 
incoming electrons. As in Fig. 3, this effect leads to a pronounced peak marking the 
magnetoconductance yielding again information about the magnetic subband structure 
at higher energies. 

For fixed kpD/Tr, it is interesting to analyse the effect of temperature and variations 
to the effective mass as exposed by various materials. To this end, we define the 
enhancement (on) - suppression (off) ratio of current fiow lon/^off in the linear response 
regime when we apply the highest magnetic field we considered, i.e., B = O.SBc- In the 
upper panel of Fig. 5, the temperature dependence of the lon/^off ratio is shown for an 
array with varying number N of coupled-dots. Remarkably enough the results hardly 
modify with > 3 in support of our previous remarks. We observe that relatively 
large ratios in excess of 100 can be achieved for temperatures up to ~ lOK and can 
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Figure 5. (Upper panel) Ratio lon/Ioff of the current flow in the on {B = 0.3-Bc ~ 
1.13T) and off {B = 0) state as a function of temperature for an array of AT = 2, 3, 5 
coupled-dots. (Lower panel) Temperature dependence of the lon/Ioff ratio for various 
effective masses m* . hrD/n and Ep are the same as in Fig. 3. 

be preserved to Ion/ 1 off > 10 for temperatures up to ~ 26K. Note that this behavior 
may be drastically improved with a selective choice of materials and geometry. A search 
in the parameter space for the latter is best done with an exhaustive analysis of the 
calculated magnetic subband structure at each field value which is beyond the scope of 
this study. Rather, in the lower panel of Fig. 5 we show that probing of the magnetic 
subband structure in materials with lighter effective mass would be greatly enhanced 
and feasible at higher temperatures. 

4. Concluding remarks 

A few additional comments regarding the effects of electronic interactions and spin are 
in order. An estimation of the ratio of the single-electron charging energy U for each 
decoupled dot over the interdot coupling t (determined mainly by D/W) yields U/t < 1 
for most dot-modes with wide energy bands forming in the array. In this case, there are 
small interaction corrections to the conductance spectrum within the first conduction 
channel (/ci^D/vr < 2) whose main effect is to sharpen the resonances observed in the 
magnetoconductance (Fig. 3), while not modifying the main transport mechanism [24]. 
At higher fcpD/n interaction effects are even less significant. Interactions may modify 
the transport profile for some distinct very narrow bands at kpn/ir < 2 for which 
U/t > 1 or if L>/VF is greatly reduced. Such an analysis of weakly coupled dots is not 
within the scope of the present study. Finally, the addition of the Zeeman splitting due 
to spin provides very small corrections to the energy and can be safely neglected for the 
magnetic fields considered here. 

To summarize, we have presented an investigation of ballistic transport through a 
finite array of coupled-dots from the perspective of a quantum mechanical magnetically 
timable mechanism that redefines bright and dark transport windows. The latter have 
been respectively identified as the energy bands and gaps of the electronic structure of 
the corresponding one-dimensional artificial crystal despite the small number of dots. 
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Thus, by tracing their magnetic field dependence we showed that the prccTirsor of 
magnetic subband formation in the energy spectrum can be readily observed. The broad 
energy range of the transport windows also reveals a well defined mechanism that yields 
magnetically controlled currents with large enhancement - suppression ratios which can 
extend up to several tens of Kelvin depending on material parameters. With present 
technology such a device can be realized within a region of ~ 300 nm at a magnetic 
field of ~ 0.5T. 
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